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Abstract 

We argue that chiral symmetry breaking in three dimensional QCD can be identified 
with Neel order in 2-dimensional quantum antiferromagnets. When operators which 
drive the chiral transition are added to these theories, we postulate that the resulting 
quantum critical behavior is in the universality class of gauged Yukawa matrix models. 
As a consequence, the chiral transition is typically of first order, although for a limited 
class of parameters it can be second order with computable critical exponents. 



One of the most intriguing features of quantum spin systems is their relationship to 
gauge theories. This connection was originally used to study chiral symmetry breaking in 
quantum chromodynamics (QCD), where the strong coupling limit resembles a spin system 
More recently, the analogy has been exploited to prove that certain gauge theories break 
chiral symmetry in the strong coupling limit [Q, ^. It has also been used to formulate 
mean field theories for magnetic systems ||^. For the most part, these works use the formal 
similarity between a gauge theory and a spin system at the lattice distance scale. Recently 
it has been suggested that the analogy is much broader in that it can account for the quasi- 
particle spectrum and other infrared features of the two systems 0. In this letter we shall 
present evidence for the latter by discussing a common feature of the phase diagrams of 
2-dimensional quantum antiferromagnets and 3-dimensional QCD. The dependence of the 
chiral symmetry breaking pattern on the number of flavors and colors of quarks in QCD is 
similar to that of the antiferromagnet where the rank of the spin algebra and the size of 
its representation play the same role as the number of flavors and colors, respectively. We 
shall also study the critical behavior associated with a chiral or Neel phase transition. Such 
a transition must be driven by operators which are added to the QCD or antiferromagnet 
hamiltonian and which have the appropriate symmetries. We argue that these transitions 
fall into a universality class which can be analyzed using the epsilon expansion. We show 
that in many cases they are fluctuation induced first order transitions. 

It has been observed that, in both 2+1-dimensional QED and QCD [^], there ex- 
ists a critical number of flavors such that if Np < Np^^' the model breaks chiral symmetry 
spontaneously and if Np > Np^^' the theory is in a chirally symmetric, deconfined phase. 
For large Np and for large number of colors Nq the equation of the critical line is approxi- 
mately Np — = 0. A heuristic argument for this behavior is that when Np » Nq 
internal gluon exchanges and the gluon self-coupling are suppressed by factors of Nc/Np. 
Resummation of leading order diagrams, which are chains of bubbles, produces an effective 
interaction which falls off like 1/r, rather than the tree level In \r\. The weak coupling of 
order Nc/Np and mild infrared behavior of this resummed theory result in a chirally sym- 
metric, de-confined phase. When Np is small, the effective coupling is large and can generate 
a condensate, which is already seen in QED [0. In fact, in QCD, when Np << Nq all pla- 
nar diagrams contribute to processes, making the effective interaction string-like |^ and the 
theory is in a confining and chiral symmetry breaking phase Numerical simulations 
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TT| of 3-dimensional QED support this scenario with A'^^'''*' ~ 4. 

Mass operators for basic 2-component fermions in 2+1 dimensions are pseudoscalars and 



break parity exphcitly |12|. Massless 2+1 dimensional QCD with an odd number of flavors 
of 2-component fermions is afflicted with the parity anomaly ||T3| which generates a parity 
violating Chern-Simons term and also fermion mass term by radiative corrections. With an 
even number of flavors, there exists a parity and gauge invariant regularization and QCD is 
the 2+1- dimensional analog of a vector-like gauge theory in 3+1 dimensions. In particular 
a kind of chiral symmetry can be deflned. It is known that, in this case, parity cannot be 



broken spontaneously and therefore to study chiral symmetry breaking it is necessary 
to seek parity conserving mass operators. Following [0, |, |15|] we shall use Np species 
of 4-component fermions. The flavor symmetry of massless QCD in this case is actually 
SU{2Np). We will add operators to the action which reduce the symmetry to SU{Np), for 
example, the gauged Nambu-Jona Lasinio (NJL) model with four-fermion interaction, 

S = Jd^x (^^F'^^ + ^^,D,^ + ^ {^T^^y^ (1) 

where is a generator of SU (Np) in the fundamental representation. The 4-fermi operator. 



which is renormalizable in the 1/Nc expansion |jT8|, can drive the chiral phase transition 
with condensate (j)^ =< ipT^ip >. The results of [|^, ^ indicate that if Np < Np^'^ the order 
persists even when A = 0. Gauged NJL models, when analyzed by solving the gap equation 



|T6| , exhibit second order behavior at a surface in the space {Np, Nc, X)- Our analysis will 



indicate that for a large range of parameters fluctuations make this transition flrst order. 
Our results do not apply to the hypothetical case where Np or Nc are varied to drive the 



transition 



The 2-dimensional generalized antiferromagnet has Hamiltonian 

-1 

^.pin = «: E E J^i^)J^iy) (2) 

<x,y> A=l 

with < x,y > nearest neighbor sites x and ?/ on a square lattice and the spin operators J^{x) 
are in an irreducible representation of the SU (Np) Lie algebra 

>(x),J^(l/)] =^r^^5(x,i/)J^(x) (3) 

When the representation at each site is a rectangular Young Tableau with m rows and A'^c 
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columns, it is convenient to represent the spin operators by the fermion bihnears 

Nc Np 

J\^) = E E ^U^)Tli^Ux) (4) 

a=l a,b=l 

The fermions have the anticommutator, 

{iljaaix),'iplf^{y)} = dabSapSix, y) (5) 
Constraints which project out the irreducible representation of the spin algebra are 

Np 

Gapix) = ^^L{x)^apix) - SafSm ~ Vx (6) 
a=l 

Qajsix) obeys the U{Nc) Lie algebra, commutes with the Hamiltonian and acts as the gen- 
erator of gauge transformations with gauge group U{Nc). 



The critical behavior of the antiferro magnet was examined by Read and Sachdev |19| 
using semiclassical methods. The only free parameters are the integers Nc and Np- Nc » 
Np is the classical limit of large representations, where the classical Neel ground state is 
stable with the staggered spin order parameter 

Nc 

l^ab = (-l)^'"^' < E '^aa{x)-^ba{x) > (7) 
a=l 

On the other hand, the limit A^^ >> A'',^ is the quantum limit where fluctuations are impor- 
tant and the system is in a spin disordered state. For both Nc and Np large, they find a 
line of second order phase transitions in the {Nc, Np) plane at Np = const. ■ A^^^ where the 
constant is a number of order one. 

The relationship between the antiferromagnet and QCD is a very close one. There is 
an argument in ref. which maps the strong coupling limit of lattice QCD onto the 
antiferromagnet with Hamiltonian (^. The lattice regularization of the QCD Hamiltonian 



uses staggered fermions pO 



H= E UMU:iMy) + h.c. + '^t^{E^{xy)y]+^^^^^ (8) 

<x,y> \ ^ A=l y □ \ □ □ / 

where the first sum is over links and the second is over plaquettes □ of the lattice. The 
gauge fields, which are unitary matrices U^y and electric field operators occupy links and 
satisfy the algebra 

E^ixy), E^{zw)] = if'^^^E^\xy)6{xy, zw) (9) 
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E^{xy),U.^,\=i5{xy,zw)U^yT^ . (10) 
The Hamiltonian is supplemented by the Gauss' law constraint 

Nf 

E E^{xy)T^^ + Y.i'\Mi^ap{x)-5^pNp/2r^Q (11) 

y&N{x) a=l 

which enforces gauge invariance. Here the first summation is over nearest neighbors of x. 

Staggered fermions have a relativistic continuum limit when their density is 1/2 of the 
maximum that is allowed by Fermi statistics and the kinetic Hamiltonian has phases which 



produce an effective U{1) magnetic flux tt per plaquette pO|, In order to obtain these 



phases in the continuum limit, we have chosen the sign of the third term in the Hamiltonian 
so that it is minimized by the conflguration of gauge flelds with the property < Hn U >= — 1. 
The constraint of half-fllling m = Np/2 is enforced by (pT]). The naive continuum limit yields 
2+1-dimensional QCD with gauge group U{Nc) and Np species of massless four component 
fermions. The full chiral symmetry only emerges in the continuum limit. On the lattice, 
staggered fermions have a discrete remnant of chiral symmetry (translation by one site) 
which forbids explicit fermion mass terms 0, A fermion mass term is a staggered 
density operator. For example, a latticization of ip{x)T^ilj[x) is obtained from the staggered 
magnetization density in eqn. (^ as Y^ab T^^Hah- Thus, the antiferromagnetic order parameter 
and the order parameter for chiral symmetry breaking with a flavor-vector condensate are 
identical. 

The argument of can be summarized as follows: The strong coupling limit, ^ oo 
suppresses fermion propagation. In the leading approximation, the Hamiltonian is minimized 
by the states which contain as little electric fleld as possible and which are compatible with 
the gauge constraint ([TT|). When Np is even [0], it is possible to solve Gauss' law with 
= 0. The occupation number of each site is Np/2 and < {—lyip'l^ipaa >= 0. This 
is a degenerate state - any gauge invariant state with Np/2 fermions has the same energy. 
Because they are required to be color singlets, this is the same set of states as occurs in 
the antiferromagnet when m = Np/2, i.e. in the representation of SU{Np) whose Young 
tableau has Nq columns and Np/2 rows. Furthermore, to resolve the degeneracy, one must 
diagonalize the matrix of perturbations. These are non-zero only at second order and the 
diagonalization problem is equivalent to solving for the ground state of the antiferromagnet 
Hamiltonian @ with k = t^/e^. Finally, since the order parameters are identical, the 
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Neel ordered states of the antiferromagnet correspond to chiral symmetry breaking states 
of QCD. Thus, the infinite couphng hmit of QCD is identical to the antiferromagnet. A 
main difference between QCD with finite couphng and the antiferromagnet is that QCD 
contains electric and gauge fields which allow a fermion kinetic energy and still retain gauge 
invariance, whereas in the antiferromagnet, the fermions are not allowed to move. One 
could regard the corrections to the strong coupling limit of QCD as the addition of degrees 
of freedom and gauge invariant perturbations in the antiferromagnet which allow fermion 
propagation. In fact, ref. p suggests even a stronger correspondence, that the additional 
degrees of freedom are generated dynamically. 

A common feature of QCD and the antiferromagnet is that, aside from Np and Nc they 
have no free parameters. We could imagine adding operators of the sort that, if their coupling 
constant is varied, it can induce the chiral transition. It is tempting to speculate that these 
transitions fall into a universality class which can take into account all such modifications, 
as long as they respect the symmetries of the theory. Here, we shall restrict our attention 
to those which lead to a Lorentz invariant continuum limit. We argue that the universality 
class is described by the 4 — e dimensional Euclidean field theory, 

S = /rf^-^^(^trV0-V0+^(|^(tr0Y + ^tr0Vit< 

+^(7 ■ V + ifi'/^ea ■ A + ifi'/^e2-f ■ TrA + I <p)^ (12) 

The scalar is an x Np traceless Hermitean matrix. The 4-component spinor i\} is an 
Np X complex matrix and is a U{Nc) gauge field. In four dimensions this model has 
Euclidean Lorentz invariance, C,P and T, discrete chiral symmetry, {ip 'y^'ip, 4> —* — 0) 
and global SU{Np) flavor. (|T2D includes all operators which are marginal when D = 4. 

The evidence that (|T2p describes the universality class comes from previous work where 
we examined a similar model where gauge couplings are absent We showed that the 

anomalous dimensions of operators computed in the model (0) with Cj = were identical 
to leading order in 1/Nc and e to those of a 2 < D < 4 dimensional four-fermi theory. That 
a 4 — e dimensional Yukawa-Higgs theory has the same universal critical behavior as lower 
dimensional four-fermi theories with the same symmetries was originally suggested by Wilson 



|22|. For the case Np = 1, where the chiral symmetry is discrete, higher order computations 



have been carried out |2^, |2^ . The results, as well as those of lattice simulations, support 
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the universality hypothesis |2^. We conjecture that (|T^) represents the universahty class 
of lower-dimensional four-fermi theories with U{Nc) gauge invariance. We shall show that, 
as a consequence, for a large range of values of (NcNp), the chiral phase transition is a 
fluctuation induced first order transition. When it is second order, critical exponents are in 
principle computable in the epsilon expansion. 

In order to analyze ([T^), we use the e expansion and the renormalization group. The 
1-loop beta functions to first order in e are. 



2 _ g 2 



e 3 iV^ - 1 2 3 2 Nl + 2NfNc - 3 
"2^ " 167r2 Nc ""'^ " 8^""'^ ^ 16iV|iVc ' 



/5e. = -2^1 ^^^^^ e, > /5e. = -2^2 + ^^e2 (13) 



We can also compute the anomalous dimensions of the scalar and the fermion field. When 
Nc and Np are large we obtain 

(27-ll)(7 + 2)2 ^ ^ 

2 (7 + 2)(27-ll)4 ^ ^ 

respectively, where 7 = Np/Nc- We shall see that the region where there can be second 
order behavior is 7 < 8.3. This result is reliable for small epsilon. 

Fixed points occur at zeros of the beta functions, I3i{g*) = 0. They are infrared (IR) stable 
or ultraviolet (UV) stable if all eigenvalues of the stability matrix, df3i/dgj\g=g*, are either 
positive or negative, respectively. The fixed points of the Higgs model (ei = 62 = y = 0) 
were analyzed by Pisarski |^ and of the Yukawa-Higgs model (ei = 62 = 0) in ref. |p!8| . 

Second order phase transitions are possible when renormalization group trajectories flow 
to an IR stable fixed point. It is at the IR stable fixed points that the conformal field theory 
with all dimensional constants except the renormalization scale are set to zero, exists. 

If there are no such fixed points, the only possible phase transition is a fluctuation induced 
first order one. Yamagishi pBf formulated a criterion for this behavior. He showed that this 
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occurs when renormalization group trajectory crosses the surface in couphng-constant space 
given by 

V,{g,y,e) = Q,i = l,2 (16) 

where Pi ((7,i/,e) = (4-e)((?i + (72)+/5i+/52 and P2 = (4-e)((?i + ^(72) +A + ^/32 depending 
on whether (72 > or (72 < respectively. When the renormahzation group trajectory 
crosses these surfaces two further conditions must be met. To ensure that the extremum 
is a local minimum, rather than a maximum, it is necessary that > for i=l,2 where 

= (4-e)(/?i + /52)+E^ + /?2 or D2 = (4-e)(/5i + ^/?2) + E. A'9/%(/?i + ^/92) 

depending on whether (72 > or (72 < 0. 

In order that this minimum has lower free energy than the trivial = 0, it is necessary 
that the couplings at that scale obey 

(7i + (72 < or (7i + NFg2/2 < (17) 

For the beta function (p!3D, there are two distinct cases. When Np < llNc/2 the four 
dimensional non-Abelian gauge coupling is asymptotically free. The solution of Pei = is 
at ci = and this fixed point is UV, rather than IR stable. In this case, there should be a 
nonperturbative behavior associated with confinement which is inaccessible to our compu- 
tation. When Np > llNc/2 the four dimensional gauge coupling is not asymptotically free. 
Pei has two zeros, the UV attractive one located at zero and the IR attractive one at 

ef = e 18 

^ 2Nf-11Nc ^ ' 

Using this solution in (3y, we see that the Yukawa coupling constant always has an UV 
attractive fixed point at zero coupling and an IR attractive fixed point at 



^ ' (llA^c - 2A^f) / + 2NfNc - 3 ^ ' 



Finally, using ([T9|), the equations for fixed points of the matrix self- couplings are 
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Numerical investigations show that IR stable fixed points (in fact only one) exist when 
(ll/2)A^c < A^F < N^iNc). The upper critical N^{Nc) intersects the line {ll/2)Nc at 
Nc = 1 and when Nq and Np are large, Np{Nc) ^ lO-TA^c-- 

Naively one would expect second order behavior for Nc and Np in this entire region. 
However, in some cases, renormalization group trajectories can satisfy the conditions for 
first order behavior before they reach the IR fixed point. In particular, if the IR fixed point 
is such that, gl + g2 < 0, then the far IR behaviour of the effective potential is approximately 
^[0] ~ {9i + g2)<P'^^^~^^^^'^~^^ and the trivial configuration = is a local maximum of the 
potential. In addition, if this {gl,g2) is the limit of a fiow originating from the ultraviolet 
stability wedge, {gi + (72 > fl gi + NFg2/2 > 0}, the potential is bounded from below, and 
since it has a local maximum at the origin, it must have a minimum when 7^ 0. Hence, in 
this case, there is always a first order phase transition, even though an IR stable fixed point 
exists. 

In fact, it is easy to see that the fiow in this case intersects the surface ([T6|). Since gl+g2 < 
0, any trajectory which fiows from the UV fixed point at the origin (in a direction within the 
stability wedge) to the IR fixed point must intersect the surface "P = in (|1^). Consider a 
trajectory which begins near the UV fixed point at the origin, where the couplings are small 
and Pi ~ —egi. Since the trajectory must be in the region + (72 > 0, P ~ (4 — 2e){gi + (72) 
which is positive. At the IR fixed point, the beta functions in V vanish and it takes the 
form P = (4 — €){gl + (72) which is negative. Since V{g, y, e) changes sign as we follow the 
trajectory from the UV to the IR fixed point, it must go through zero at least once. 

We have verified numerically for a {Nc, Np) = {(2, 15), (2, 18), (100, 1000)}, with e = 0.1 
that they all satisfy the criterea for a first order phase transition. Generally, IR fixed points 
with gl + g2 < occur in the range Np{Nc) < Np < Np{Nc)- In this region, the theory 
has an IR stable fixed point, but undergoes a first order phase transition. For large Np and 
Nc, Np{Nc) = 8.3Nc. 

The other possibility is that the IR stable fixed point has the property gl + g2 > 0. In this 
case, we have found that there are two kinds of trajectories, those which fiow from the UV to 
IR fixed points without encountering the surface V = and lead to second order transitions, 
and those which visit the region where gi + g2 < 0. The latter trajectories cross the surface 
twice and satisfy the conditions for first order behavior at one of the intersections. In fig.l 
we have plotted several RG fiows for the case A'",^ = 2, Np = 13 and e = .1. (The features 



8 



of the flow are insensitive to e in the domain .1 < e < 1. Near e = 1, the relevant coupling 
constants are large and one would not expect the 1-loop approximation that we have used 
to be accurate there.) The stars indicate that the flow has crossed the surface P = in the 



region defined by > and (|T7|). Hence the transition associated with those trajectories 
is first order. The crosses denote points where the flow intersects the surface (|T^), but not 
when Qi + g2 < 0. They do not correspond to global minima of the effective potential and 
therefore do not indicate a first order transition. 

The physical quantum spin j antiferromagnet corresponds to Np = 2 and Nq = 2j. 
With A''^; = 1 (j = 1/2) the non-abelian field must be removed, and the resulting theory 
has no IR fixed point, indicating a first order transition. Np = 2, Nc > 2 (i.e. j > 1) is in 
the asymptotically free regime and hence these antiferromagnets cannot be analyzed by our 
techniques. We speculate that confinement is associated with a nonperturbative IR fixed 
point of the gauge coupling. In that case, it is likely that these antiferromagnets would have 
a second order transition. 
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Figure caption: Fig. l:Projection of the renormalization group trajectories in five dimen- 
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RG Flow 

Nc = 2Nf = 13eps = 0.1 




